W4 thermochemistry of P2 and P4. Is the CODATA heat of 
formation of phosphorus atom correct?* 

Amir Karton and Jan M. L. Martin* 
Department of Organic Chemistry, Weizmann Institute of Science, 
IL-76100 Rehovot, Israel. Email: comartin@weizmann.ac.il 
(Dated: Received May 24, 2007; Accepted June 28, 2007.) 
( Mol. Phys. (Peter Pulay festschrift) TMPH-2007-0140) 

Abstract 

The high-accuracy W4 computational thermochemistry protocol, and several post-W4 methods, 
have been applied to the P2 and P4 molecules. Contrary to previous studies, we find the 
experimental thermochemistry to be fundamentally sound. The reaction enthalpy for P4 — > 
2P2 has a very significant contribution from post-CCSD(T) correlation effects. We derive a 
gas-phase heat of formation for the phosphorus atom of AH^ [P(g)]=75.54±0.1 kcal/mol and 
Aii^ 298 [P(g)]=75.74±0.1 kcal/mol, in the upper half of the CODATA uncertainty interval. 



* Dedicated to Professor Peter Pulay on the occasion of his 65th birthday. 
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I. INTRODUCTION 



The most important component of phosphorus vapor is tetrahedral P4. It is also a 
building block of the most stable solid allotropes (white and red phosphorus), as well as an 
important component of the liquid phase 1]. 

P4 is a challenging molecule to describe well in terms of both 1-particle basis set and 
n-particle space ("electron correlation treatment"). Its peculiar bonding pattern has been 
studied in some detail by Ahlrichs et al.[2]. Perhaps the most comprehensive ab initio 
study to date has been the work of Persson, Taylor, and Lee (PTL) 3j, who studied the 
thermochemistry, geometry, and anharmonic force field of the molecule at the MP2 (second- 
order many-body perturbation theory) and CCSD(T) (coupled cluster with all singles and 
doubles plus quasiperturbative triples [4] ) levels with basis sets of up to spdfg quality Their 
work buttressed the earlier contention of Haser and Treutler^] that the established bond 
distance 6| is considerably too long due to misassignment of the rotational fine structure of 
the Raman spectrum. Subsequent IR work[?J confirmed this conclusion. 

In addition, the PTL anharmonic force field calculations showed implausibly large 
discrepancies with the experimental band origins, which were resolved by subsequent 
remeasurements 8[] . 

Finally, PTL addressed the endothermicity of the reaction P 4 — > 2P 2 , and found a 
significant discrepancy with experiment, a substantial part of it resulting from a very 
large core- valence contribution. Haworth and Bacskay (HB)]9|, in a later benchmark 
thermochemical study on a number of phosphorus compounds including P2 and P4, were 
able to rationalize the core- valence contribution as an artifact of the basis set used. However, 
despite their work being carried out near the CCSD(T) basis set limit and including 
both core- valence and scalar relativistic corrections, HB's best calculated total atomization 
energies (TAE values) still imply a reaction energy for P 4 — > 2P 2 that is almost 2 kcal/mol 
more endothermic than experiment [101] . 

This is not a purely academic issue. This selfsame reaction represents the weakest link in 
the determination of the CODATA 11[ gas-phase heat of formation of the phosphorus atom. 
(The other two links, namely the heat of vaporization of P4 and the dissociation energy of 



P2, are fairly precisely known 



11 



12j.) Said heat of formation enters the equation whenever 
one calculates the heat of formation of any phosphorus compound — be it by ab initio, 
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semiempirical, or density functional techniques. 

Very recently, we developed a next-generation, fully ab initio, computational 
thermochemistry method known as W4 theor v[13l ] , th at was shown to be able to reproduce 



the best ATcT (active thermochemical tables 14|, [15j, ll6|]) atomization energies for several 
dozen small molecules with an RMS deviation of 0.08 kcal/mol (95% confidence interval 
of 0.16 kcal/mol). Among other things, this method was applied effectively in a recent 



revision 



171 | of the atomic heats of formation of boron, aluminum, and silicon. It would 



appear to be an appropriate tool for trying to resolve the P4 issue. In the present paper, we 
will show that the experimental thermochemistry in general, and the CODATA AHj [P(g)] 
value in particular, are in fact fundamentally correct. 

II. COMPUTATIONAL METHODS 



All calculations reported in the present work were performed on the Linux cluster of the 
Martin group. All SCF, CCSD and CCSD(T) calculations were carried out using MOLPRO 



2006.l[21[. Some CCSDT calculations were carried out using the Austin-Mainz-Budapest 

22j program system; the remaining post-CCSD(T) calculations were 



version of the ACES II 

carried out using an OpenMP parallel version of Kallay's general coupled cluster code 



MRCC 



23 



24 



251 ] interfaced to the Austin-Mainz-Budapest version of ACES 1122]. SCF 



and correlated DBOCs were obtained using the relevant module 26j of the PSI3 open source 
quantum chemistry code 27]. All basis sets employed belong to the correlation consistent 



family of Dunning and coworkers 



28 
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30 



311 ] . Convergence of iterative processes was 



accelerated by means of Pulay's direct inversion of the iterative subspace (DIIS) method 32]. 

The computational protocols for W4 theory, for the lower-cost methods W2.2, W3.2, and 
W41ite, and for the post- W4 methods W4.2, W4.3, and W4.4 have been documented in great 



detail elsewhere 
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18j |. and will not be repeated here. For the sake of clarity in the following 



discussion, we do note that the W4 energy consists of the following components (summarized 
in Table I): SCF, valence CCSD correlation, valence parenthetical triples (T), higher- 
order connected triples T 3 — (T), quasiperturbative connected quadruples (Q), higher-order 
connected quadruples T 4 — (Q), connected quintuples, CCSD(T) level inner-shell correlation, 
scalar relativistic effects (second-order Douglas-Kroll-Hess approximation 19]), first-order 
spin-orbit coupling, and diagonal Born-Oppenheimer correction (DBOC). W4.2 and higher 



3 



add a contribution for core- valence higher-order triples, while W4.3 uses more extended basis 
sets in calculating and extrapolating the post-CCSD(T) correlation effects (besides adding 
a connected sextuples term). W4.4 theory |l_8| employs improved extrapolations, as well as 
core-valence parenthetical quadruples and a correlation correction to the DBOC. 

III. RESULTS AND DISCUSSION 
A. At the CCSD(T) limit 

All relevant thermochemical data can be found in Table I. Convergence of various 
contributions as a function of basis set is summarized in Table II. 

The valence-correlation-only CCSD(T)/cc-pV(Q+d)Z bond distance for P 4 , 2.1987 A, 
is nearly 0.01 A longer than the all-electron CCSD(T)/cc-pwCVQZ bond distance, 2.1893 
A. Combining this with the computed ro — r e of PTL, 0.005 A, this implies r =2.194 A, 
in excellent agreement with the experimental value of Boudon et al. 7j], ro=2.1958 A. (An 
earlier measurement pj of 2.2228±0.0005 A was previously rejected by PTL on the basis of 
their calculations, which were subsequently confirmed by the Boudon et al. work.) 

The valence-only and all-electron calculated bond distances of P2 are 1.9019 and 1.8941 
A, respectively. The latter value is in excellent agreement with the experimental one from 
Huber and Herzbergp^, r e =1.8934 A. 

It was recently noted, in the context of a basis set convergence study on higher-order 
—on e ffeCtS Q, th- the use of valence-only correlated reference geometries has 
insignificant consequences for first-row systems, but has more noticeable effects for second 
row compounds with strong bonds. In the present work, we find a small but non-negligible 
0.05 kcal/mol for P2 and a whopping 0.13 kcal/mol for P4 — the latter is not surprising in 
light of the hefty geometry change. 

The zero-point vibrational energy of P4, 3.91 kcal/mol, was derived from combining the 
anharmonicity constants of PTL with the revised experimental fundamentals of Ref. f|. 
The one for P 2 , 1.11 kcal/mol, derives from the spectroscopic constants given in Huber and 
Her zb erg [3]. 

At the W4 level, the valence CCSD(T) contributions to TAEo add up to 116.09 kcal/mol 
for P2, and 289.32 kcal/mol for P4. These values are lower by 0.6-0.7 and 1.1-1.5 
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kcal/mol, respectively, than the corresponding limits of HB. (The ranges given refer to the 
spread among the various basis set combinations extrapolated from.) Of these respective 
discrepancies, 0.12 (P2) and 0.24 kcal/mol (P4) reflect the different definitions of the 
restricted open-shell CCSD(T) energy used for phosphorus atom: Watts- Gauss-Bartlett [33], 
a.k.a. RHF-UCCSD(T), in the present work, and Knowles-Hampel- Werner 341 ] . a.k.a., RHF- 
RCCSD(T), in HB. The remainder may somewhat reflect the different reference geometries 
used (HB employed B3LYP/6-31G(2df) density functional calculations): stretching or 
compressing bonds causes a redistribution of binding energy between valence and core- 
valence contributions that is much more significant than the overall change in the sum of 
these parts. A more important component of the difference probably results from the fact 
that while HB extrapolated the total valence CCSD(T) energy, the present authors carried 
out separate extrapolations of the SCF, valence CCSD correlation, and (T) contributions, 
and in addition extrapolate singlet-and triplet-coupled pair correlation energies separately 
as advocated by Klopper[35[. Furthermore, we use more extended aug-cc-pV(5+d)Z and 
aug-cc-pV(6+d)Z basis sets. Using smaller basis sets with our procedure — as is done in 
W2.2 and W3.2 theory — changes the atomization energies of P2 and P4 by just -0.07 and 
+0.08 kcal/mol, respectively, suggesting that our W4-level extrapolation should be close to 
Nature. 

Likewise, our core-valence contributions differ somewhat from those of HB: for P 4 we 
obtain 1.76 kcal/mol extrapolated from CCSD(T)/aug-cc-pwCVTZ and aug-cc-pwCVQZ 
basis sets, while their result with the rather small cc-pCVTZ basis set is 2.2 kcal/mol without 
counterpoise correction, and 1.4 kcal/mol with counterpoise correction. HB's core- valence 
contributions for P 2 likewise bracket ours. 

Our scalar relativistic contributions, -0.28 (P2) and -0.72 (P 4 ) kcal/mol, are somewhat 
larger in absolute value than those of HB (-0.2 and -0.5 kcal/mol, respectively): while 
they employ first-order Darwin and mass- velocity corrections [36] from a CASPT2I37I wave 
function with the G31arge basis set 38], we use the second-order Douglas-Kroll-Hess[l^] 



approximation from a CCSD(T) wave function with a relativistically contracted aug-cc- 
pVQZ basis set {39]. The latter approximation should normally be the more rigorous one. 
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B. Beyond the CCSD(T) limit 

We shall now consider post-CCSD(T) contributions. The %[(T)] diagnostic for P 2 , 
8.4%, indicates a system with moderate, bordering on severe, nondynamical correlation. 
(The largest T2 amplitude in the CCSD calculation, 0.10 per degenerate component, is 
related to the 7r — > ir* excitation. For additional diagnostics, see Table III.) As a result, error 
compensation between higher-order triples (-1.00 kcal/mol at the W4 level) and connected 
higher-order excitations (+1.46 kcal/mol connected quadruples, +0.10 kcal/mol connected 
quintuples) is less than perfect, and any CCSD(T)-based thermochemistry protocol would 
significantly underestimate -Do- I n fact, the difference between W4 and the higher- level 
W4.3 protocol, 0.20 kcal/mol, is the largest we have thus far encountered. The difference 
breaks down as 0.06 kcal/mol higher-order triples, 0.06 kcal/mol connected quadruples, 
0.03 kcal/mol connected quintuples, 0.03 kcal/mol from the newly introduced connected 
sextuples term, and 0.03 kcal/mol from higher-order triples in the inner shell correlation 
term. Applying our highest-level protocol, W4.4 theory, only causes an additional 0.02 
kcal/mol change in Dq: decreases in the CCSD and (T) extrapolated contributions are 
compensated by a 0.11 kcal/mol contribution from core- valence connected quadruples. 

Our W4.4 A), 116.22+0.10 kcal/mol at the CCSD(T)/cc-pV(Q+d)Z geometry and 
116.27+0.10 kcal/mol at the CCSD(T)/cc-pwCVQZ geometry, at first sight seems 
significantly higher than the spectroscopic value of 116.06+0.09 kcal/mol. The latter (from 



a 1932 Herzberg paper[40j]) derives from a predissociation limit in the C X S+ state to 
P ( 4 S) +P ( 2 D) , where it is not clear whether the latter is in the J = 3/2 or J = 5/2 
fine structure level. Taking the fine structure data from the NIST database, we obtain 
116.10+0.07 kcal/mol in the former, and 116.06+0.07 in the latter case. Especially the 
J = 3/2 case agrees to within overlapping uncertainties with our calculations. 

We considered yet larger basis sets for calculating the post-CCSD(T) corrections, namely 
cc-pV5Z for T 3 — (T) and (Q), cc-pVQZ for T 4 — (Q), and cc-pVTZ for T 5 . Results can be 
found in Table II. The higher-order triples contribution would change by -0.01 kcal/mol 
relative to W4.3 but by +0.01 kcal/mol relative to W4.4. The combined quadruples 



contributions would decrease by 0.04 kcal/mol, but in our experience [18j, this basis set 
expansion tends to be mostly cancelled by the effect of expanding the basis set for the 
core-valence higher-order triples — which, incidentally, is the reason (besides computational 
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intractability) why all three expansions are neglected in W4.4 theory. Unfortunately, a 
CCSDT/cc-pwCVQZ calculation is not feasible (yet) for a molecule with two second-row 
atoms. Finally, the quintuples contribution would go up by 0.03 kcal/mol. 

Let us now turn to P4. The highest-level calculation we were able to perform was W4 
theory, and even that stretched our computational resources (8-core AMD Opteron and Intel 
Cloverton machines with 16 to 32 GB of RAM and 2 TB of striped disk space each) to the 
limit. With a %[(T)] diagnostic of 7.5%, the CCSD(T) method is significantly in error here 
too; yet, we note that no single T 2 amplitude is larger than 0.05 for P 4 . We find a whopping 
-3.28 kcal/mol from higher-order triples, 2.57 kcal/mol from connected quadruples, and 0.13 
kcal/mol from connected quintuples. The computed TAEo of 285.96±0.16 kcal/mol at the 
standard W4 reference geometry, and 286.09±0.16 kcal/mol at the core-valence reference 
geometry, agree well with the Gurvich TAE of 285. 9±1 kcal/mol, but this is a rather 
hollow victory in light of the stated uncertainty. We expect that the computed TAE would 
go up significantly if we were able to perform W4.3 or W4.4 calculations. 

Interestingly, for P4 CCSD(T) errs on the opposite side from the case of Dq(P2), such 
that the error on the reaction energy P4 — > 2P 2 actually gets amplified: overall, CCSD(T) 
overestimates the full CI limit by 1.70 kcal/mol for this reaction. 

The experimental value actually derives from D (P2) and the rather uncertain measured 
enthalpy for P 4 — > 2P 2 , evaluated in Gurvich et al. [lOj as 53.83±1 kcal/mol. As can be 
seen in Table I, our computed reaction energy for P4 — > 2P 2 converges fairly rapidly with 
the level of theory, and it can reasonably be assumed that the W4 value, 53.97 kcal/mol at 
the standard reference geometry and 54.00 kcal/mol at the core- valence geometry, is not far 
from Nature. Combining the latter with our W4.4 D (P 2 ) would suggest that W4 and W4.4 
TAE values would differ by about 0.4 kcal/mol, which is not impossible but does seem 
excessive. 286.53±0.10 kcal/mol would seem to be a plausible upper limit, and 286.10±0.16 
a plausible lower limit. Intermediate values of 286.12 and 286.20 kcal/mol are obtained by 
combining the computed reaction energy with the two possible experimental _Do(P 2 ) values. 
Our best estimate, 286.32±0.22 kcal/mol, splits the difference between our estimated upper 
and lower limits. 

Combining one-quarter of our best estimated TAE [P 4 (g)] with one-quarter of the quite 



precisely known [10], 1 1 II ] Aif£ [P 4 (g)]=15.83±0.07 kcal/mol, we obtain an estimated gas- 
phase heat of formation of phosphorus atom, Aif^ [P(g)]=75.54±0.1 kcal/mol. Using the 
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CODATA heat content functions of P(g) and white phosphorus (see footnotes to Table 
I), we obtain AiJ^ 298 [P(g)]=75.74±0.1 kcal/mol, which is in excellent agreement with the 
CODATA reference value of 75.65±0.24 kcal/mol, but carries a much smaller uncertainty. 
Our calculations do seem to suggest that the true value is in the upper half of the CODATA 
uncertainty interval. 



IV. CONCLUSIONS 



Summing up, the high-accuracy W4 computational thermochemistry protocol, and 
several post-W4 methods, have been applied to the P2 and P4 molecules. Contrary to 
previous studies, we find the experimental thermochemistry to be fundamentally sound. The 
quite significant contribution of post-CCSD(T) correlation effects to the reaction enthalpy 
for P4 2P2 illustrates the importance of their proper treatment in accurate computational 
thermochemistry work. We derive a gas-phase heat of formation for the phosphorus atom 
of 75.54±0.1 kcal/mol, in the upper half of the CODATA uncertainty interval. 
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TABLE I: Component breakdown of total atomization energies of P2 and P4, as well as reverse 
dimerization reaction. Computed and CODATA heats of formation at absolute zero of phosphorus 
atom. All data in kcal/mol 





P 2 - 
2 P 


P 4 
4 P 


P 4 - 
2P 2 


Aff? [P(g)] 


SCF 




W4 anH W4 y 

VV t: dllvA VV rt.-A. 


38.93 


121.40 


43.54 




CCSD 


aug-cc-pV({5,6}+d)Z b 


W4,W4.2,W4.3 


67.45 


146.17 


11.27 






anp'-rr-nVf'-f 5 6ll-d'lZ c 


W4.4 


67.43 


146.12 


11.27 




\ 1 ) 


micr rr rA/f-fO ^1-4-^^7^ 
dug p v qVi'Jj 1 


W4 W4 9 W4 3 


9.71 


21.82 


2.40 








W4.4 


9.67 


21.70 


2.36 






rr-nV-fD T\7J 


W4 W4 2 W4 3 

VV Tt) VV l.i^ VV rt.U 


-1.00 


-3.28 


-1.29 






rr-nV-fT QlZ-f 


W4.3 


-0.94 










cc-dVAT Q1Z» 


W4.4 


-0.96 








T, 


1 1 xrr-nVTZ,' 1 

l.X AW> U V _1_ Zj 


W4 W4 2 

V V -±; V V rfc. j. 


1.46 


2.57 


-0.34 






pp nV-TT 017* 


W4 3 W4 A 

VV 'i.O) VV rt-rt 


1.52 








^ 


cc-pVDZ(no d) 


W4 W4 2 

VV .1 - VV rr-ii 


0.10 


0.13 J 


-0.07 






pp rWT>7 


W4 3 W4 4 

VV 'i.O) VV rt-rt 


0.13 








Tr 
1 f) 


pp nVD7fnn rH 


W4 3 W4 A 

VV rt.ij) VV l.rt 


0.03 










rTSTVn /antr rr nw^V-fT Cl\7, 


W4 


0.73 


1.76 


0.31 




shell 


_A \_ J v_y O -L' X / v^L, t_/ W v_y V J- ZJ 


W4.2 3 


0.76 










-lAOl /rr nwflVTZ 


W4.4 


0.87 








i'A 1 'ifn,- 


DK rCSTVT'l /ancr rr nVOZ 

J_ / IV v_, V_j O J_ / I _1_ J / alig LL. p V ZJ 


W4 and W4 y 

VV txIH_l VV ^ . A 


-0.28 


-0.72 


-0.17 




DBOC fe 


WTT/niio- pp r>VT7 


W4 Wi 9 WzL ^ 
VV 4:, VV Z, vv^±.o 


0.01 


0.02 


0.00 






4-CTSD /rp-nVDZ 


W4.4 


0.01 








A-M l 




W4 only 


0.02 


0.04 


0.00 




ZPVE 




W4 


1.11 


3.91 


-1.68 




AE n 




W2.2 


115.36 


286.60 


55.87 


75.61 m 






W3.2 


115.41 


285.11 


54.28 


75.23 






W41ite 


115.48 


285.03 


54.07 


75.22 






W4 


116.00 


285.96 


53.97 


75.45 






W4.2 


116.02 












W4.3 


116.20 






[75.55]™ 






W4.4 


116.22 






[75.56] n 


Ageom. 




all 


0.05 


0.13 


0.03 


0.03 






best 


116.27±0.10 


[286.32±0.22]° 




[75.54±0.1] n 


Experiment 




CODATA 


116.0±0.7 


285.9±1.0 


53.87±0.7 


75.45i0.24 1 ? 






Gurvich 


116.06±0.09 p 




53.83±1 


75.45±0.249 



The notation cc-pV{X,Y}Z, for instance, indicates extrapolation from cc-pVXZ and cc-pVYZ basis sets. 
a extrapolated using the Karton-Martin^S] modification of Jensen's extrapolation formulaj^ 

b T\ term from largest basis set, singlet- and triplet-coupled pair energies are extrapolated using the two-point A + B/L° 
expression of Halkier et al.[45j with a =3 and 5, respectively. 

c same expressions, but using Schwcnke's extrapolation which is equivalent to a =3.06967 for singlet pairs and 
a =4.62528 for triplet pairs. 

d extrapolated using the two-point A + B/L 3 expression of Halkier et al.j^l 



(Footnotes continued on next page of manuscript.) 
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e same expression, but using Schwenke's extrapolation 46], which is equivalent to a =3.22788 

f extrapolated from AVQZ and AV5Z basis sets using the two-point A + B/L 3 expression of Halkier et al. [45! 

9 extrapolated using empirical A + B/L 5 / 2 found in Ref. |is| 

h Approximated as 1.1 X (Q)/cc-pVTZ + 1.1 X [Q-(Q)]/cc-pVDZ, as proposed in Ref. [13] 

1 (Q) contribution extrapolated cc-pV{T,Q}Z, Q-(Q) contribution obtained as difference between CCSDTQ and CCSDT(Q) 

calculations with the cc-pVTZ basis set. 

•» Value for P 4 is approximated by CCSDTQ5 A -CCSDTQ 

k Post-Hartree-Fock corrections to the DBOC as obtained at the CISD/cc-pVDZ level are found to be +0.003 kcal/mol for 
P2 and +0.001 kcal/mol for P4. This correction exhibits very weak basis set sensitivity: A. Karton, B. Ruscic, and J. M. L. 
Martin, J. Mol. Struct. (THEOCHEMj 811, 345-353 (2007). 

1 differences between ACES II and MOLPRO definitions of the valence CCSD(T), extrapolated from PVTZ and PVQZ basis 
sets using the two-point A + B/L 3 expression of Halkier et al.j^l 

m Computed AHj Q [P(g)] values are derived as (l/4)(AflJ Q [P 4 (g)]+TAE [P4(g)]), where the experimental 

AH° [P 4 (g)] =15.83+0.07 kcal/mol is obtained from CODATA AH° 2g8 [P 4 (g)]=58.9±0.3 kj/mol, 

[#298 - -ffo][P4(g)]=14.10±0.20 kj/mol, and [H 2 98 - H ] [P(cr,white)]=5.360±0.015 kj/mol. Note that using the best 

currently available fundamental frequencies ;8j for P4, the rigid rotor-harmonic oscillator [H29S — Ho] [P4(g)]=14.046 kj/mol. 

71 Assuming that Ai? r ,o for the dimerization reaction, 2 P2 — >P4, is converged at the W4 level. Further basis set effects then 

amount to (l/2)Afl o(Pa)- 

° See the text. 

p Spectroscopic value extracted from Herzberg[4o|, Value derives from predissociation of C x Sj state to P( 4 S3/2)+P( 2 -D x), 
where x G {3/2, 5/2}. Combining Hcrzberg's predissociation limit of 51969+24 cm -1 with the latest atomic excitation energy 
data from the NIST database[41| ( 2 D 3 / 2 : 11361.02 cm' 1 and 2 D 5/2 11376.63 cm" 1 ), we obtain 116.10+0.07 kcal/mol if 
2 D 3/2 and 116.06+0.07 kcal/mol if 2 D 5/2 - 

1 From CODATA AHj 298 [P(g)]=316.5±l kj/mol, [H 2 gs - H ] [P(g)]=6. 197+0.001 kj/mol, and 
[H 2 98 - i?o][P(cr,white)]=5.360±0.015 kj/mol. 
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TABLE II: Basis set convergence of various contributions to the total atomization energy (kcal/ 







P 2 








P 4 








AV(T+d)Z AV(Q+d)Z AV(5+d)Z AV(6+d)Z 


AV(T+d)Z AV(Q+d)Z AV(5+d)Z AV(6+d)Z 


SCF 




JO. / z 


QQ QQ 
JO.OO 


^8 Q1 


120 24 




121 17 


121 36 


CCSD T-pair 


23.86 


24.45 


24.58 


24.61 


66.64 


68.94 


69.50 


69.66 


CCSD S-pair 


36.07 


40.27 


41.73 


42.32 


59.95 


69.79 


73.54 


75.01 


CCSD Ti 


-0.29 


-0.30 


-0.31 


-0.31 


-0.58 


-0.60 


-0.62 


-0.62 


CCSD total 


59.65 


64.42 


66.01 


66.62 


126.01 


138.13 


142.43 


144.05 


(T) 


9.09 


9.36 


9.53 


9.59 


19.79 


20.82 


21.31 






PVDZ 


PVTZ 


PVQZ 


PV5Z 


PVDZ 


PVTZ 


PVQZ 


PV5Z 




-0.07 


-0.72 


-0.85 


-0.90 


-0.78 


-2.54 






(Q) 


1.04 


1.43 


1.57 


1.61 


1.77 


2.66 






2i-(Q) 


-0.11 


-0.15 


-0.17 




-0.32 








T 5 


0.13 a 


0.16 






(a) 










ACVTZ 


ACVQZ 






ACVTZ 


ACVQZ 






core- valence 


0.67 


0.70 






2.00 


1.86 







(a) With PVDZ(no d) basis set: 0.10 kcal/mol for P 2 , 0.13 kcal/mol for P 4 . 



TABLE III: Diagnostics for importance of nondynamical correlation 





%TAE[SCF] a 


%TAE[(T)] a 


%TAE 


%TAEp4 4 


-T 5 ] a 


Ti[47j 


D i[48] Largest T 2 NO occupations 








[post-CCSD(T)] a 






diag 


nostic amplitudes HDOMO 6 LUMO 
















- CCSD(T)/cc-pVTZ - 


p 2 


33.2 


8.29 


0.48 


1.33 




0.018 


0.033 0.118 (x2) 1.899 0.085 


P4 


41.9 


7.53 


-0.20 


0.93 




0.018 


0.037 0.048 1.916 0.061 



"Percentages of the total atomization energy relate to nonrelativistic, clamped-nuclei values with inner shell electrons 
constrained to be doubly occupied, (from W4 theory) 
b Highest doubly occupied molecular orbital. 
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